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Collisional interaction limits between dark 
matters and baryons in 'cooling flow' clusters 
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ABSTRACT 

Presuming weak collisional interactions to exchange the kinetic energy between dark 
matter and baryonic matter in a galaxy cluster, we re-examine the effectiveness of 
this process in several 'cooling flow' galaxy clusters using available X-ray observa- 
tions and infer an upper limit on the heavy dark matter particle (DMP)— proton 
cross section cTxp. With a relative collisional velocity V^— dependent power-law form 
of (Jxp = CTo(T^/10'^ km s"^)" where a < 0, our inferred upper limit is UQ/m^ ^ 2 x 

10~^^ cm^ GeV""'^ with m^ being the DMP mass. Based on a simple stability analysis of 
the thermal energy balance equation, we argue that the mechanism of DMP— baryon 
collisional interactions is unlikely to be a stable nongravitational heating source of 
intracluster medium (ICM) in inner core regions of 'cooling flow' galaxy clusters. 

Key words: cooling flows — cosmology: theory — dark matter — galaxy: clusters: 
general — radiation mechanisms: general — X— rays: galaxies: clusters 
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1 INTRODUCTION 

Astrophysical and cosmological measurements together with 
numerical simulation experiments indicate that the cold 
dark matter (CDM) constitutes most of the matter in the 
Universe, even though the fundamental physical nature of 
such dark matter particles (DMPs) remains unknown. In 
the simplest scenario, these cold DMPs are presumed to 
be collisionless and they interact with each other or with 
other baryons only through the mutual gravity at the present 
epoch. While the collisionless CDM model is successful in ex- 
plaining the formation of large-scale structures in the Uni- 
verse, observational contradictions with numerical simula- 
tions appear inevitable regarding structures on sub-cluster 
scales, e.g., the prediction of a higher number of dwarf galax- 
ies than that actually observed. Spergel & Steinhardt (2000) 
revitalized the concept of strongly interacting massive par- 
ticles (SIMPs) to confront this issue and suggested a self- 
interacting cross section per unit DMP mass a^^/ruy^ ~ 
10~^* cm^ GeV~^ where Uxx is the collisional cross section 
among DMPs and mx is the DMP mass. Along this line of 
reasoning, if DMPs are strongly self-interacting, then simi- 
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lar strong interactions would also be equally expected to ex- 
ist between DMPs and baryons. Elastic scatterings between 
DMPs with nuclei could generate a recoil energy of the order 
of ~ 10 kcV, which would then be detectable in underground 
and underwater particle experiments. Such 'direct detection' 
experiments currently limit the cross section to the order of 
10"*^ — 10"'"^ cm^ for a cold DMP of mass in the range of 
~ 10- 10^ GeV (e.g., Akerib et al. 2004). Limits on such cob 
lisional interactions in a wider mass range were investigated 
by various physical experiments and observations of several 
astrophysical processes, such as /3/3 decays, cosmic-ray de- 
tections, the galactic-halo stability, the cooling of molecu- 
lar clouds, proton decay experiments, the existence of old 
neutron stars and the Earth (e.g., Starkman et al. 1990), 
satellite experiments in space (e.g., Wandelt et al. 2001), 
primordial nucleosynthesis and cosmic rays (e.g., Cyburt et 
al. 2002), and cosmic microwave background anisotropy and 
large-scale structure power spectrum (e.g., Chen et al. 2002). 
These experiments and observations have provided a limit 
on the mass-dependent cross section per unit DMP mass be- 
tween DMPs and protons as Uxp/m-x < 10^^** — lO^^** cm^ 
GeV"^ for a DMP mass mx > 1 GeV. 

Recent high-resolution observations of the hot X-ray 
emitting gaseous intracluster medium (ICM) by Chandra 
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(e.g., Peterson et al. 2001) and XMM-Newton (e.g., Kaas- 
tra et al. 2001; Tamura et al. 2001) satellites have revealed 
deficits of cool gases (with gas temperatures much less than 
the virial temperature Tvir) in the core of so-called 'cool- 
ing flow' galaxy clusters, inconsistent with predictions of the 
conventional radiative cooling models (e.g., Cowie & Binney 
1977; Fabian & Nulsen 1977; Mathews & Bregman 1978; 
Stewart et al. 1984; Fabian 1994 and extensive references 
therein). Several heating mechanisms have been proposed 
to resolve this 'cooling flow dilemma', such as the inward 
thermal conduction from hot outer regions (e.g., Narayan & 
Medvedev 2001 and references therein), energy injections as- 
sociated with central activities by an active galactic nucleus 
(e.g., Churazov et al. 2002), and outward or inward acoustic 
wave heating (e.g., Pringle 1989; Fabian et al. 2003a, b; Fu- 
jita et al. 2004; Feng, Zhang, Lou & Li 2004). The resonant 
excitation of internal gravity modes ((;— modes) in the ICM 
by orbiting galaxies was explored by Balbus & Soker (1990); 
they examined the processes of excitation, propagation, am- 
plification, damping of such galaxy cluster g—modes in the 
context of providing thermal energy in 'cooling flow' galaxy 
clusters. 

Qin & Wu (2001) proposed coUisional interactions be- 
tween baryons and heavy DMPs {m^ ^ rup with rrip be- 
ing the proton mass) as a major nongravitational heating 
mechanism for ICM in the core of a galaxy cluster. As- 
suming DMPs and baryons have comparable velocity dis- 
persions, the kinetic energy of a single DMP would then 
be much larger than that of a baryon because of the pre- 
sumed mass difference m^ 3> m-p. Therefore in numerous 
elastic collisions, kinetic energies of DMPs can be systemat- 
ically transferred to the baryonic ICM to balance the radia- 
tive cooling in X-ray bands by hot electrons. By equating 
heating and cooling rates, they estimated a specific cross 
section of a^p/m^ ~ 10~^^ cm^ GeV~^ for m^/mp > 10^. 
This is actually a requirement of compensating the radia- 
tion cooling and should not be regarded as a limit of any 
sort. It should also be noted that they took strong lensing 
cluster CL0024-I-1654 to infer the DMP- baryon interaction 
cross section for balancing the cooling, which may not be 
proper for the following reason. The significant discrepancy 
by a factor of 3 ~ 4 between the mass profiles derived from 
X-ray observations and gravitational lensing effects shows 
that the ICM may not be in a static equilibrium and might 
be still collapsing (e.g., Kneib et al. 2003; Zhang et al. 2005), 
and the bimodal velocity distribution of cluster galaxies in- 
dicates a merger of two systems with a mass ratio of 1 to 2 
(e.g., Czoske et al. 2002). The main heating mechanism of 
CL0024-f 1654 could well be the gravitational collapse. 

Chuzhoy & Nusser (2006) re-considered the ICM heat- 
ing scenario of Qin & Wu (2001), corrected their calcula- 
tions and derived a similar cross section for the heavy DMP- 
proton elastic coUisional interaction. They found that, if CTxp 
is independent of the relative velocity V of colliding parti- 
cles, a thermal equilibrium state between heating by DMPs 
and radiative cooling by hot electrons of ICM would be al- 
ways unstable. However, in galaxy clusters with T > 2 keV, 
a stable energetic balance may be achieved for a relative 
velocity F— dependent cross section a^p ex V^ with a < —3. 

There are two major simplifications in both analyses 
of Qin & Wu (2001) and Chuzhoy & Nusser (2006). First, 
they took the temperature and density of the ICM and the 



density of DMPs for typical values, rather than the relevant 
distributions determined by X-ray observations with high 
angular and spectral resolutions. Secondly, they estimated 
the velocity dispersion of DMPs either simply similar to that 
of baryons or by the results of numerical simulations, rather 
than a dynamically self-consistent 'true' value obtained by 
solving the Jeans equation (Binney & Tremaine 1987; Sub- 
ramanian 2000; see also Ikebe et al. 2004 for the case of 
galaxy cluster A1795). Since the specific cross section limit 
they derived (at the centre or at the virial radius) is highly 
sensitive to the chosen parameters [see equation (8) of Qin 
& Wu (2001)], it is crucial to investigate this ICM heating 
mechanism more carefully using an actual sample of 'cooling 
flow' galaxy clusters, well observed in X-ray bands by the 
XMM-Newton and Chandra satellites in space. 

This paper is structured as follows. In §2, we set the 
upper limits for DMP-proton elastic coUisional cross section 
per unit DMP mass. In §3, we demonstrate that the DMP- 
proton coUisional interaction alone is unlikely to be a stable 
ICM heating mechanism to compensate the radiative cool- 
ing of ICM. Discussion and conclusions are contained in §4. 
Details on the X-ray cooling function can be found in Ap- 
pendix A. 

In our theoretical model consideration, we have adopted 
the currently favoured standard ACDM cosmology with the 
cosmological parameters h = 0.7, Sim = 0.3 and JIa ~ 0.7 in 
conventional notations. 



2 DARK MATTER-BARYON COLLISIONAL 
INTERACTION CROSS SECTIONS 

2.1 Model Description 

The X-ray radiative cooling rate of the ICM can be inferred 
observationally and may be approximately represented by 
the cooling function A (i.e., radiative energy loss rate per 
unit volume) in terms of temperature T, baryon mass den- 
sity pb and abundances Z of the ICM as described in Ap- 
pendix A. The radiative 'cooling time' t^ of a galaxy cluster 
is then defined by (e.g., Sarazin 1988) 
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where e is the ICM internal energy per unit volume, fee is 
the Boltzmann constant, fi is the mean molecular weight of 
the ICM, and irip is the proton mass. The cooling time tc is a 
function of radius r in general. For 'cooling flow' clusters, the 
estimated central cooling time is shorter than their cosmic 
age, or the Hubble time tn = Hg with Ho being the Hubble 
constant. The so-called 'cooling radius' re is defined as the 
radius such that tc{rc) — tn- 

The thermal conduction in ICM across magnetic field 
may be negligible (e.g., Cowie & Binney 1977; Stewart et 
al. 1984; Sarazin 1988; Fabian 1994) due to smaU gyro- 
radii given by rg = jmv±c/{ZceB) where 7 = [1 — {v?, + 
v'j_)/c^]~^'^ is the relativistic Lorentz factor, m is the par- 
ticle mass, W|| and v± are the particle velocity components 
parallel and perpendicular to the local magnetic field B, Z^e. 
is the particle electric charge, and c is the speed of light r] 

^ By numerical simulation analysis, Tao (1993) argued however 
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For thermal electrons and an average |-B| ~ 1/iG in a typical 
ICM, we estimate Vg to be a few thousand kilometers. The 
central regions of 'cooling flow' clusters will cool down sub- 
stantially within a timescale of to since their formation; in 
the absence of other effective heating mechanisms to com- 
pensate the radiative loss in X-ray bands, the hot ICM core 
would then collapse under self-gravity. Although the ICM 
behaviour in 'cooling flow' clusters is inhomogeneous on 
smaller scales due to thermal instabilities (e.g., Field 1965; 
Mathews & Bregman 1978; Malagoli et al. 1990), their main 
large-scale properties may be grossly modelled by smooth 
subsonic flows in a theoretical description (e.g., Fabian et al. 
19841"] Assuming a quasi-spherical symmetry and a quasi- 
steady state, the ICM evolution may be described by the 
following two equations, namely, the mass conservation of 
baryons 



M = 4:TTp\jVr 



(2) 



and the energy conservation of baryons in ICM (see Ap- 
pendix B) 



PhV 



K + (j) + 



= A-H 



(3) 



where M is the baryon mass accretion rate, the ICM ther- 
mal pressure p = phk^T / (fimp) follows the ideal gas law, 
(j} is the total gravitational potential (including that of the 
dark matter halo), v'^/2 is the baryon kinetic energy per 
unit mass, and H represents the possible heating function 
(i.e., the heating rate per unit baryon volume). The specific 
baryon enthalpy /ib in energy conservation equation (|3| is 
given by 
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(7 - l)pb 2fimp 



(4) 



where the ratio of specific heats 7 is taken to be 7 = 5/3 
for the ICM. The estimated typical 'cooling flow' speed v ~ 
10 km s"'^ is much less than the ICM sound speed and thus 
1)^/2 is ignorable as compared to hh in energy equation (|3|. 
We could have added a magnetic term {_B^)/(47rpb) 
within the parentheses on the left-hand side of energy equa- 
tion ([3| with {B'i) being the mean square of a completely 
random magnetic field Bt transverse to the radial direc- 
tion (Yu & Lou 2005; Yu et al. 2006; Wang & Lou 2007; 
Lou & Wang 2007). For a diffuse ICM random magnetic 
field (Bt)^''^ of strength ~ 10~^G away from the core re- 
gion (e.g., Clarke et al. 2001; Carilh & Taylor 2002) and 



a typical proton number density of rip 



10^ 



(e.j 



Sarazin 1988; Voit 2005), we estimate a typical Alfven speed 
of ~ 100 km s~^ much less than the ICM sound speed. Such 
a magnetic field strength may give rise to an anisotropic 
distribution of electrons but may not be significant in the 
sense of bulk flow dynamics. By the equipartition argument, 
a random magnetic fleld (i?^)^'^ rnay reach strengths as 
strong as ~ 10 — 30 /xG in the ICM core region of some 



that 'tangled' magnetic field may not be effective enough to sup- 
press the thermal conduction in haloes of galaxy clusters. 
^ Chuzhoy & Nusser (2006) adopted different equations to de- 
scribe the ICM behaviour. In particular, they ignored ICM flows 
as well as the effect of significant gravitational heating. In these 
two aspects, our model appears more general and realistic. 



galaxy clusters (e.g., Fabian 1994; Hu & Lou 2004). De- 
pending on the actual proton number density in the range 
of Up ~ 10~^ — 10~^ cm~^, magnetic pressure is less than 
(e.g., Dolag & Schindler 2000) or may be comparable to the 
thermal pressure in the core region for relaxed clusters of 
galaxies. For our main purpose of inferring the upper limit 
of DMP mass ttIx based on data of galaxy clusters, it suffices 
to examine the purely hydrodynamic case. 

Now we consider the possible ICM heating mechanism 
due to collisions between DMPs and baryons. As usual, we 
assume Maxwellian velocity distributions for both protons 
and DMPs with (one dimensional) velocity dispersions Vp 
and Ux respectively. For the case of V independent cross 
section, we have recalculatecrj the energy transfer rate per 
unit volume through elastic collisions as 



H = 8 



1 + 



1/2 



1 



m^n^n^v^ax 



ipVp/{m„v^) 
(1 + mp/m^y 



(5) 

where rip and rix are the number densities of protons and 
DMPs, respectively. For DMPs much heavier than protons 
(i.e., nix 3> wip) and as v^ ~ rip in a gravitationally virialized 
system, the last factor in the form of a division on the right- 
hand side of equation S approaches unity. 

For the case of V dependent cross section o^p, we take 
the case that the cross section has the form of a power law, 
namely a^p = ao{V /VoY . To avoid complicated calcula- 
tions and as an example of illustration, we only consider the 
case of heavy DMPs. Adopting expression (9) of Chuzhoy & 
Nusser (2006), we would then have 

H^Qx 2''mpnpnxi;^+"Vo-Vo(l + r^p/r'x)^'"^"'''' (6) 

as an approximate heating function. 

In our model anaylsis, the scattering cross section cTxp 
between a DMP and a proton is a free parameter to be con- 
strained by several available X-ray observations of 'cooling 
flow' galaxy clusters. The scattering cross section between 
DMPs and helium nuclei UxHc is expected to be 4crxp for 
incoherent scatterings, or 16(Txp for coherent scatterings, or 
for the spin-dependent case (e.g., Chen et al. 2002). For 
simplicity, we shall take CTxHc = 0, noting that the other two 
alternatives would change the results only slightly. 

The mass density and temperature distributions of the 
ICM can be inferred from X-ray observations independently. 
The one-dimensional velocity dispersion of ICM protons is 
Vp = {k^T/mpY''^ . Under the approximations of quasi- 
spherical symmetry, quasi-hydrostatic equilibrium and ideal 
gas law, the total enclosed cluster mass distribution obeys 
the following condition 



GMr 



knT 
firup 



dlnT dlnpb 



dr 



dr 



(7) 



where Mr is the total enclosed cluster mass (dark matter 
and baryonic matter together) inside radius r, G = 6.67 x 
10~* g~^ cm^ s~'^ is the gravitational constant. The mass 
density of DMPs can then be inferred from 

1 dAIr 



Px(r-) = 



inr^ dr 



Pb(r) . 



(8) 



•^ We have verified that the constant coeflicient in equation (4) 
of Qin & Wu (2001) was in error. 
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In most model calculations for clusters of galaxies, the 
DMP mass density distribution is fitted with the universal 
NFW mass profile obtained by numerical simulations (e.g., 
Navarro, Frenk & White 1995, 1996, 1997), namely 



Px(r-) = Pxo('-/rs) ^(1 + r/vsY 



(9) 



where rs is a radial scale and pxo is a DMP mass density 
scale when r = 0.48 r^. For r <^ Ts, DMP mass density p^{r) 
scales as r~^, while for r 3> Ts, px{r) scales as r~^. There 
are also other possible dark matter mass density profiles 
Px(r) ex r-'{r + r,)'"" with 1 < / < 1.5 and 2.5 < q < 3 for 
different combinations of inner and outer radial scalings of 
DMP mass density (e.g., Moore et al. 1998, 1999; Jing 2000; 
Rasia et al. 2004; Voit 2005). 

In a spherical hydrostatic equilibrium with no mean 
streaming motions such that Vr = ve = ^4, = 0, the ve- 
locity dispersion of DMPs in a galaxy cluster and the mass 
distribution is related by the Jeans equation (e.g., Binney & 
Tremaine 1987; Subramanian 2000) 

— (px«,,) + — Px«. +Px^3-=0. (10) 

Here, the velocity anisotropy parameter Pa is defined as 



Pa = l-v^(r)/[2vf.{r)] , 



(11) 



where Vt = Vg+ v^ and Vr, ve, v^ are the radial and angular 
velocity dispersions with respect to the mean velocity which 
is almost zero in a quasi-hydrostatic equilibrium. The formal 



solution of Jeans equation ( 10 1 with a constant Pa is 



2f \ "^ 



-■213a 



P^{r) 



„2/3, 



X p^{x)GMxAx 



(12) 



Numerical simulations indicate a variation of the anisotropy 
parameter Pa between /3a = at the halo centre and 
Pa ~ 0.6 at the virial radius (e.g., Colin et al. 2000). For 
Pa as a function of r, Jeans equation ( 10 1 can be integrated 
numerically with additional approximations (e.g., Subrama- 
nian 2000). As the form of Pa cannot be determined directly 
by observations, we simply set Pa = for an isotropic ve- 
locity dispersion in our model analysis of galaxy clusters. In 
fact, our calculations indicate that a nonzero Pa only affects 
v^ at a level of < 10%. 



2.2 The upper limit of elastic collisional cross 
section constrained by X-ray observations 

Observations have shown that the ICM temperature is com- 
parable to the equivalent temperature of cluster galaxies 
(e.g., Jones & Forman 1984), implying that the ICM has 
not changed its temperature very much from the initial hy- 
drostatic equilibrium state since the epoch of cluster for- 
mation. We thus consider a more general problem of ICM 
heating and cooling [see equations (pi and (|3|] in a quasi- 
hydrostatic manner. 

Integrating equation (|3| from the centre r = to the 
cooling radius Vc and using equation H, we derive 



h + (t>+Y 



A-H 



M/(47rr2) 



dr 



(13) 



For 'cooling flow' galaxy clusters, we have 0o < ^('"c) and 
To < T{rc) according to satellite X-ray observations. Near 
the cluster centre, w ~ 10 km s~^ is very much lower than 



the sound speed there. If we take either a quasi-hydrostatic 
or a quasi-magnetostatic solution (Lou & Wang 2006, 2007) 
for cluster evolution just before a complete virialization or 
equilibrium, then u ~ at r = and thus < v^{r^)/2. 
By the above consideration and estimates of three terms at 
both r = Tc and r = separ ately , we show the left-hand 
side (LHS) of energy equation (13 1 is positivePjIn fact, this 
requirement may be relaxed a bit. In fact, as long as the 



LHS of energy equation ( 13 1 is non- negative and irrespective 
of the relative magnitudes of each individual terms at both 
r = re and r = 0, we have for the baryonic ICM the following 
inequality for the net heating rate 



(A - H)r^dr > , 



(14) 



and the collisional cross section is limited from above by 

a^p< f " Kr^Ar I f " Hir^dr , (15) 

where Hi = H/a^p according to expression ([5|, or by 



O-qVo " < 



HiV^dr 



(16) 



where H2 = HVo/ao according to expression m. By the 
above analysis and reasoning for the roles of flow kinetic 
energy and magnetic energy and in reference to inequali- 



ties ( 14 1 — ( 16 1, our model calculations and X-ray data com- 
parisons can be much simplified within a purely quasi- 
hydrostatic framework and that is what we do in the fol- 
lowing. 

We calculate the DMP-proton elastic collisional cross 
section limit using X-ray observations of a sample of five 
presumably relaxed galaxy clusters with apparent 'cooling' 
cores, namely, Abell 478 (e.g., Pointecouteau et al. 2004), 
AbeU 1795 (e.g., Ikebe et al. 2004), AbeU 1983 (e.g., Pratt 
& Arnaud 2003), AbeU 1991 (e.g., Pratt & Arnaud 2005), 
and PKS 0745-191 (e.g., Chen et al. 2003; Pointecouteau 
et al. 2005). Physical parameters estimated for this galaxy 
cluster sample are summarized in Table 1. This sample of 
galaxy clusters covers wide ranges of redshift 2 (0 — 0.1), 
cluster virial mass M200 ([1 - H] x lO" Mq), and ICM tem- 
perature T (2 — 7 keV), respectively. These chosen galaxy 
clusters are analysed with the deprojection technique based 
on the X-ray data taken from the XMM-Newton EPIC to 
reveal the real spectra of gaseous ICM in different spherical 
shells and to determine the deprojected temperature and 
the mass distribution of the gas in a galaxy cluster. The 
XMM-Newton EPIC is the most sensitive X-ray space tele- 
scope which also has high spatial and spectral resolutions, 
and thus meets all the requirements for a detailed spectral 
analysis. The statistical study on the comparison of cluster 
mass measurements using strong gravitational lensing and 
using X-ray observations shows an excellent mutual agree- 
ment for 'cooling fiow' galaxy clusters, which are perceived 
as dynamically more relaxed systems and the quasi-spherical 
hydrostatic assumption used in X-ray mass determinations 
should be valid (Allen 1998; Wu 2000). On this ground, the 
dark matter mass density distribution used in this paper is 

"* Note that in reference to a large radius, h + if)-\-v'^ /2 should be 
negative for a gravitationally bound system, while the difference 

{h + <j) + v^ jl^r^ — {h + (/> + ^^/2)o can be positive as shown. 
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Table 1. Estimated physical parameters for a sample of five 'cooling flow' clusters of galaxies 



Cluster 


2 


T (keV) 


C200 


-R200 (kpc) 


Afaoo (IO^^Mq) 


re (kpc) 


M (M© yr- 


ija 


A478 


0.0881 


6.8 


4.22 ±0.39 


2060 ±110 


10.8 ±1.8 


17-1 + 113 
-'-'-'--115 


73g+ii4 

'•50„434 




A1795 


0.0616 


5.7 


4.47 ±0.27 


1760 ± 30 


6.54 ±0.35 


129+1° 


■^^^-213 




A1983 


0.0442 


2.3 


3.83 ±0.71 


1100 ± 140 


1.59 ±0.61 


34+11 


g„+10.8 
■^•^-6.0 




A1991 


0.0586 


2.6 


5.78 ±0.35 


1106 ±41 


1.63 ±0.18 


52+_% 


37+?6 




PKS0745 


0.1028 


7.0 


7.05 ±0.28 


1880 ± 130 


8.34 ±0.84 


i26j;^^ 


^70+399 





"^ The dark matter mass deposition rates from 'cooling flows' (column 8) are taken from White et al. (1997). 
While their model calculations used cosmological parameters of h = 0.5 and go = 0-5> the corresponding 
parameter values here are all comparable to theirs. 
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Figure 1. On the basis of data analysis and for the case of ex- 
pression ([5I, the upper limits of speciflc cross section crxp/nJx for 
the flve selected 'cooling flow' galaxy clusters (namely, PKS 0745- 
191, AbeU 478, Abell 1795, Abell 1991, and AbeU 1983) are shown 
here. The value range of the estimated upper limit for Cxp/rrax is 
mainly due to uncertainties of the cooling radius Tc and the virial 



Figure 2. For the relative velocity V— dependent case of expres- 
sion ijml, the upper limits of cross section o-o/rrix as functions 
of exponent a for the five selected 'cooling flow' galaxy clusters 
(namely, PKS 0745-191, AbeU 478, Abell 1795, Abell 1991, and 
Abell 1983) are shown here by flve separate solid curves with 
corresponding galaxy cluster names labelled explicitly. 



regarded as reliable, especially for the lensing galaxy cluster 
PKS 0745-191 (e.g., Allen 1996). 

For the velocity independent case as given by expression 
([5|, the upper limits of the specific cross section a^^/m^ for 
DMP-proton collisions are displayed in Figure 1. The limits 
of cTxp/jTix for the five selected sample galaxy clusters lie in 
the range of [2.0, 4.9]xl0~^^ cm^ GeV~^ and a common 
consistent value would be taken as (Txp/mx < 2 x 10~^^ 
cm^ GeV~^. For the velocity dependent case as given by 
expression (|6|, the value of cro/mx is displayed in Figure 2. 
Here, we take the mean value of cooling radius Tc of every 
cluster for simplicity. The upper limit of DMP specific cross 
section is cTxp/rrix < 2 x 10"^'^(F/10^ km s"^)" cm^ GeV"^ 

Note that (Txp/mx is constant only for m^ ^ m^ in 
expressions (|5| and (Im. Our estimated upper limit applies 
well to rrix > 10'' GeV. The cross-section limit for the lighter 
DMPs (rrix = 10 — 10^ GeV) has been strictly constrained 
by direct detection experiments (e.g., Wandelt et al. 2001). 



3 DARK MATTER PROTON COLLISIONS AS 
THE ICM HEATING SOURCE IN 'COOLING 
FLOW CLUSTERS OF GALAXIES 

The energy equilibrium equation of the ICM can be written 



dln5/di = dlne/dt = {H - A)/e 



(17) 



2/3n 



where S = kBT/(iimpp^ ) is the ICM entropy (e.g., Voit 
2005). We note that the ratio Vp/v^ in two expressions ([5| 
and m is almost a constant in dynamically and thermally 
well relaxed galaxy clusters (Hu & Lou 2007). Hence for a 
specific 'cooling flow' cluster, H oc pbT" and A ex p^T^ , and 



equation (17 1 becomes 



A\n{T/pl'^)/dt = AT'^ - Bp^T^ 



(18) 



where A and B are two coefficients independent of ph and 
T. A systematic comparison between the ICM temperatures 
T and velocity dispersion of DMP base on the X-ray obser- 
vations indicates that T could only slightly change (AT <1 
keV) in the cosmic age of clusters (Hu & Lou 2007), which 
is consistent with the numerical simulations of the cluster 
evolution (e.g., Bryan & Norman 1998). If the RHS of equa- 
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tion ( 18 1 maintains zero, i.e. cooling is fully compensated 
by heating, we have dlogpb ~ — MlogT, where |6| < 1 for 
T ~ lO'^ — 10* K. Therefore pb will be roughly a constant in 
such a case. 

Now we fix the value of T and perform an isothermal 



stability analysis on the solution pb of equation (18 1. Sup- 



pose the DMP-proton elastic coUisional heating can com- 
pensate radiative cooling for certain values T = To and 
Pb = Po, let us consider the stability of pb in equation (171 
at a constant temperature. A small baryon density variation 
5pb obeys 



2 dSpb 
3pb dt 



BT'-^Sp^ 



(19) 



A simple estimation indicates that Sph will double in a 
timescale of 



2"" X loVr 



103 km s-i , 

Px o-o/mx 



(20) 



,0.02 M© pc-3 10-25 cm2 GeV-i , 
Therefore, without other heating sources, the solution of 



equation (171 appears unstable for whatever values of a, at 



least in the case of a fixed ICM temperature profile. 



4 DISCUSSION AND CONCLUSIONS 

Based on two different elastic collision models and X-ray 
observations of galaxy clusters, we have re-examined effects 
of coUisional interactions between heavy DMPs and protons 
in 'cooling flow' clusters of galaxies and estimated a more 
reliable upper limit for the specific cross section between 
DMPs and protons by solving the equation of energy conser- 
vation for five "cooling flow" galaxy clusters using available 
X-ray data. In the regime of m,x S> rrip, the upper limit for 
c"xp ~ foC^/lO^ km s"^)" with a < can be expressed as 
o"o/?7ix < 2 X 10~^^ cm^ GeV~^, which is fully consistent 
with the earlier results. Similar to other astrophysical con- 
straints, this upper limit is independent of the underlying 
model for particle physics governing coUisional interactions. 

In our model calculations, we have assumed that the 
total gravitational potential (j) remains invariable during the 
process of galaxy cluster evolution. The gravitational po- 
tential of a galaxy cluster is dominantly determined by 
the amount of DMPs. Besides the weak energy-momentum 
transfer between the DMPs and baryons, the density distri- 
bution of dark matter is only affected by the mild accumu- 
lation of cool gas in the central region of clusters. The exact 
solution of such complex co-evolution of DMPs and baryons 
should be treated with Fokker-Planck equation (e.g., Binney 
& Tremaine 1987). We will pursue the relevant problems in 
contexts of galaxy clusters in separate papers. 

We have also demonstrated an intrinsic instability in 
explaining the problem of 'cooling flow' galaxy clusters via 
the heat transfer by heavy DMP-proton coUisional scatter- 
ing. According to expression ([5| , the DMP-proton coUisional 



heating of ICM works only when mpVp < 



In 'cooling 



flow' galaxy clusters, Vp/v^^ ~ p ~ 0.6 (Hu & Lou 2007) 
and thus m^ > 0.6 GeV is required in this model. However, 
as the coUisional cross section in the DMP mass range 0.5 



GeV< mx < 10^ GeV is less than 5 x 10"^** cm^ (e.g., Wan- 
delt et al. 2001), the light DMP coUisional heating mecha- 
nism can then be completely ignored. Therefore, we rule out 
the possibility of DMP coUisional scattering as a major sta- 
ble non-gravitational heating sources of the ICM in 'cooling 
flow' galaxy clusters. 
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APPENDIX A: X-RAY COOLING FUNCTION 

The radiative energy loss from a hot plasma per unit vol- 
ume and per unit time is a cooling function of temperature, 
density, and chemical abundances of various elements. The 
cooling functions for a hot plasma under equilibrium con- 
ditions have been calculated over a photon energy range of 
0.001 — 30 keV and a range of abundances for 0—1.0 Zq 
by Sutherland & Dopita (1993). They defined the normal- 
ized cooling function by Ao = A/(ncni), where ric and Ui are 
the electron and ion number densities, respectively. Tozzi & 



Table Al. Parameters for the Normalized Cooling Function Ao 



Z{Zq) 


Ci 


Ca 


Cs 


ric/n.H 


Ui/Uh 


M 





-0.003 


0.0605 


0.0204 


1.128 


1.064 


0.58 


0.1 


0.0193 


0.0632 


0.0218 


1.131 


1.064 


0.58 


0.32 


0.480 


0.0658 


0.0306 


1.165 


1.080 


0.60 


1.0 


0.1434 


0.0762 


0.0355 


1.209 


1.099 


0.62 



Norman (2001) used an approximate analytic expression to 
fit Ao in the following polynomial form 

Ao = [CiiksT)-'--'' + C2{k^T)°-'^ + C3] X 10"'' erg cm^ s"' , 

(Al) 
where k^T is in unit of keV. The three constant coefficients 
Ci (i = 1,2,3) and n^jn^, ni/nn, [i, depend on the metal- 
licity Z. We have refitted these parameters and summarized 
the results in Table \KA\ for the convenience of reference. For 
an arbitrary metallicity Z, we take the parameters by the 
following polynomial interpolations: 

Ci = -0.003 



C2 -- 


= 0.0605 + 0.034i 


C3 -- 


= 0.019 -H0.041Z 


ric/nu = 


= 1.12 + 0.14Z- 


m/riH = 


= 1.06 -f o.oez - 


M = 


= 0.575 -f- 0.074Z 



0.26Z-0.41Z^-h0.29Z^ , 
0.069^^ -f 0.052Z^ 



- 0.025Z^ 
0.05Z^ , 
0.02Z^ , 

- O.OSOZ'^ 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 



Expression (All with parameters defined by equations 



(A2) — (A7) can grossly reproduce the X-ray cooling function 
of Sutherland & Dopita (1993) to within a few percent in 
the typical ICM thermal energy range of 1.0 keV < k^T < 
10 keV. 



APPENDIX B: DERIVATION OF ENERGY 
CONSERVATION FOR ICM BARYONS 

The energy conservation equation for ICM baryons (equa- 
tion 3) is derived in this Appendix B. We start with the 
thermodynamic relation 

d/ib = Pb 'dp + dg = Pb Mp + Pb '(^ - A)di , (Bl) 

where /ib is the enthalpy per unit baryon mass, Q is the 
heating rate per unit baryon mass, p is the ICM (baryon) 
gas pressure, and A and H are cooling and heating rates per 
unit volume, respectively. Differentiating along the moving 



direction of a bulk baryon flow, equation ( Bl I can be written 
as 

v-(pbVp) = v-Vftb+Pb'(A-i/) . (B2) 

Using the Euler equation (e.g., Landau & Lifshitz 1959) 



9v 

at 



+ (v ■ V)v = -p^^Vp - V<f) 



(B3) 



where is the total gravitational potential (including that 
of a dark matter halo), and a vector identity 

..2 \ 

V X (V X v) , (B4) 



(v-V)v = v( y 



equation ( B2 1 appears as 



9v 
v.-.-vV 



+ 4> + hh] + 



(H-A) 



Pb 



(B5) 
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Here we assume that the moving direction of a cooling flow 
is radially inward (i.e., negative), such that 

-^ = ^[^+<t> + hb]+ ^ '- ■ B6 

at or \ 2 J phv 

We now estimate and compare the magnitudes of terms in 



equation (B6l. For example, 





dv 


V 


10 km s~^ 


. 










dt ~ 


t^ 


10"' yr 


3 X 10 m 


s 




(Hi) 


d(t) 


GM, 


- > 


GMr 


10" MqG 


' 10" 


"m 


s-2 


dr 


r^ 




r^ r^iJaoo 


(1 Mpc)2 






(B8) 


d 


(-) 


r^ 


v'^ (lOkm/s)^ 


- ~3 X 10~" 


m s' 


-2 


(B9) 


dr 


V^/ 




2r IMpc 










dhb 
dr 


r 


5 (IkeV) 

- r^ 

2ump IMpc 


~ 10~"m s 


-2 




(BIO) 



where GMr/r^ is a decreasing function with increasing r for 
NFW mass profiles obtained through numerical simulations. 
The value of {H — A)/{pf,v) term is to be determined, yet it 
should be mu ch larger than terms d{v'^ /2)/dr and dv/dt in 



equation (B6l. Therefore the LHS and d{v'^ /2)/dr term on 
the RHS of equation ( |B6[ ) may be neglected and we finally 
derive equation (3) as 



Pb«^( y +<^ + /ib) =A-//. (Bll) 



